For: "Mathematics of Materials: Percolation and Composites", Editors: 
K.M. Golden, G.R. Grimmett, R.D. James, G.W. Milton, and P.N. Sen. 
The IMA Volumes in Mathematics and its Applications (Springer- Verlag, 1997). 
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Abstract. Scaling limits of critical percolation models show major differences be- 
tween low and high dimensional models. The article discusses the formulation of the 
continuum limit for the former case. A mathematical framework is proposed for the di- 
rect description of the limiting continuum theory. The resulting structure is expected to 
exhibit strict conformal invariance, and facilitate the mathematical discussion of ques- 
tions related to universality of critical behavior, conformal invariance, and some relations 
with a number of field theories. 
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1. Introduction. Incipient percolation clusters have attracted atten- 
tion as objects of interesting physical and mathematical properties, and 
potential for applications. An example of a setup in which they play a 
role is an array of conducting elements, placed at random in an insulating 
medium, with the density adjusted to be close to the percolation thresh- 
old. In such arrays the current is channelled through fractal-like sets. The 
concentration of the current, or stress/strain in other similar setups, may 
result in high amplification of non- linear effects. The phenomenon is of 
technological interest, and plays a role in high-contrast composite materials 
and non- linear composites, utilized in thermistors and other devices [1,2,3]. 
Studies of the relevant random geometry have yielded interesting geometric 
concepts such as the celebrated (but often misunderstood) Incipient Infi- 
nite Cluster (IIC). The topic was reviewed from a physics perspective in 
an article (Stanley [4]) which appeared in Volume 8 of this series, in the 
proceedings of a workshop held at IMA in 1986. 

It is somewhat surprising that percolation threshold phenomena are 
still a source of delightful and new observations, since the subject seemed 
to be reaching its maturity already ten years ago. Nevertheless, the subject 
has recently enjoyed renewed attention; in part because it was realized that 
some entrenched notions need correction (ref. [5-13]), and in part because 
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it was realized that the scaling (continuum) limit has interesting proper- 
ties, e.g., conformal invariance, and its construction presents an interesting 
mathematical challenge [14-16, 10]. 

This article focuses on issues related to the scaling limit of the Incipient 
Spanning Clusters (ISC). We discuss a mathematical framework for the 
direct description of the limiting continuum theory, which is applicable 
to models below the upper-critical dimension. Within it, we encounter 
some interesting fractal structures and questions related to: universality of 
the critical behavior, conformal invariance and relations with certain field 
theories. The purpose is to describe some recent developments and point 
directions for possibly interesting progress. 

2. The incipient infinite cluster and the incipient spanning 
clusters. A good starting point for the discussion of the scaling limit is 
the conceptual difference between two related terms: the Incipient Infinite 
Cluster (IIC), and the Incipient Spanning Clusters (ISC) (interpreted here 
as in ref. [12]). 

The Incipient Infinite Cluster (IIC) is a thought-provoking notion 
which has often been used in the discussion of different aspects of criti- 
cal percolation phenomena ([4]). It has provided a useful and stimulating 
concept, but one whose different uses have led to some confusion and misin- 
terpretation of theoretical predictions, in particular concerning the question 
of uniqueness (see the Stauffer paradox discussed in [12]). Consequently, 
it was recently proposed to differentiate between the two related notions 
mentioned above: the IIC and the ISC. The distinction is most apparent in 
the limits in which the two are described by different mathematical entities. 

The Incipient Infinite Cluster is a random infinite cluster (or, more 
completely, a random environment exhibiting an infinite cluster), con- 
structed by a limiting process which provides the "local" picture of the 
large clusters seen at the percolation threshold — viewed from the per- 
spective of one of their sites. This notion covers three alternative paths 
which have been explored towards the construction of the IIC: 

i. Condition on the origin being connected distance L away, and let 
L — > oo. 

ii. Raise the percolation density above the percolation threshold, con- 
dition on the origin belonging to an in infinite cluster, and then let 
P \ Pc (where the percolation density presumably vanishes). 

iii. Generate the probability distribution for a random cluster (or a 
random environment) by centering a typical random configuration 
relative to one of the sites on its spanning clusters [or just large 
clusters] . The sampling is to be done with equal weights over all the 
sites connected to the boundary of [0, L] d (and then let L — ► oo), 
possibly with a corrective exclusion of a boundary zone. 

The first two procedures were pointed out by Kestcn, who proved their 
convergence and equivalence in d = 2 dimensions [17]; the third is more 
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convenient for numerical studies, and was used in [18]. Either way, one 
may see that the mathematical construction covers the microscopic view 
of the incipient infinite cluster from the view of one of the rare sites which 
at p = p c are connected very far - on the microscopic scale. (A somewhat 
different conception of the IIC is presented in ref [19]. Presumably all will 
agree in the further scaling/continuum limit). 

The Incipient Spanning Clusters, on the other hand, are naturally 
viewed on the "macroscopic" scale. They are simply the large clusters 
which reach across the finite sample (in the above example), and connect 
opposite boundary segments, as indicated in Figure 1. In order to see the 
entire collection of spanning clusters, we need to keep track of events on 
the scale of the sample, the relevant limit being: lattice spacing (a) — > 0. 




Fig. 1. Incipient Spanning Clusters - a schematic depiction of the the macroscopic 
view. 

A particular difference which stands out, and one which has caused 
extensive confusion and discussions, is that the the IIC typically shows a 
single dominant (infinite) cluster, whereas the ISC quite generally exhibit 
multiple clusters of comparable size, [12]. (A more complete discussion of 
the related Stauffer paradox is found in ref [12]). 

The continuum limit enables a natural discussion of the enhanced sym- 
metry. The highest symmetry emerges at the critical point, for which there 
is strong evidence of conformal invariance. Considerations related to 2D 
conformal fields have led to proposals for differential equations which de- 
termine some of the properties of the critical measure ([15,16]). Thus, the 
continuum limit of the spanning clusters may remind one of the Brown- 
ian motion: an object arising from physics, with fascinating mathematical 
properties, high degree of symmetry, and relation to interesting differential 
equations. 

While we focus here on ICS, similar considerations can be applied to 
the entire collection of the connected clusters which are visible on the 



4 



MICHAEL AIZENMAN 



macroscopic scale. The entire ensemble's scaling limit has been termed 
the percolation web [12]. 

3. The microscopic view; three convenient models. There is a 
variety of situations in which conducting elements are placed at random 
in an insulating medium. These elements may conduct electric current, or 
may serve as passages for a liquid seeping through a solid. The individual 
resistors, or cracks, are visible on the Microscopic Scale. Their relative 
density is our (dimcnsionless) control parameter. Following are three con- 
venient models, which offer different advantages as possible starting points 
for the construction of a (common?) continuum limit. 




Fig. 2. The microscopic view; different possibilities for the short-distance structure. 

Bond-percolation. The model is formulated on the lattice aZ d . The 
randomness is associated with the bonds (pairs of neighboring lattice sites) , 
which are occupied at random, independently with probability < p < 1 
(the control parameter). The occupied bonds are regarded as connecting. 
They may also represent conductors of conductance a > 0. 

When a bond is not occupied, its dual obstructing object can be re- 
garded as realized. That leads to the self-duality of bond percolation on 
Z 2 , which is helpful ([20,21]). In three dimensions, the dual model is that 
of random plaquettes, which may form encapsulating surfaces obstructing 
the connections ([22]). 

Droplet percolation. This models is formulated over the continuum. The 
conducting regions consist of randomly distributed d-dimcnsional balls of 
radius a (Poisson process, with possible overlaps), with density pd d x for 
the centers of the conducting cells. The relevant dimensionless control 
parameter is proportional to the density of the conducting regime: 



(3.1) 



p = p a d 
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The 2D bond model is self dual, while the droplet model is rotation 
invariant. The following model exhibits both features (a fact noticed in- 
dependently also by Benjamini and Schramm ([23]), and it, therefore, is 
our favored starting point for the construction of the purportedly common 
scaling limit. 

Voronoi-tessellation percolation. Starting from a randomly generated 
configuration of points in R d , described by a Poisson process with density 
pd d x, the plane is partitioned into the cells of the corresponding Voronoi 
tessellation. The cells are conducting, or not, independently with probabil- 
ities {p, 1 — p}. Alternatively stated, two discrete random sets (A, B C R d ) 
are generated with Poisson densities ppd d x and (1 — p)pd d x 1 and the con- 
ducting regime consists of those sites of the continuum which are closer to 
A than to B. The short-distance scale in this model, a, is related to p as in 
cq. 3.1 (it is of the order of the mean diameter of the Voronoi cells). In two 
dimensions the model is self-dual, and the critical value for p is p c = 1/2. 

An generalization we shall mention in Section 10 consists of models 
with a density profile of the form p t {x)d d x with pt(x) = t ■ g(x), g(x) 
continuous and positive. Letting t — > oo we find a scaling limit for which 
the density profile shows persistent inhomogeneity on the macroscopic scale. 
However, since the inhomogeneity corresponds to just different rates of 
approach towards a common limit, we expect it to have no visible effect on 
the continuum limit considered here. 

Planarity is a very helpful property even without strict self duality. 
Using it one can prove more for 2D models, at or near p c ([21,24]), than 
what is known about other dimensions 2 < d < 6, and in particular about 
d = 3. 

4. The macroscopic perspective. The focus of our discussion is on 
the geometric features which are visible on the Macroscopic Scale in a 

systems whose short scale structure is any of the above. Correspondingly, 
we chose the scale for our discussion so that the object occupies a fixed 
continuum-scale region Ai = [— 1, l] d , or more generally A^ = [— £,£} x 
[— 1, l] d_1 , and we let the short distance scale be a << 1, eventually taken 
to 0. 

When the sample is placed between two conducting plates, with differ- 
ent electric potentials applied to the two opposite faces: 

(4.1) <9A_ = {x e dA : Xl = -£} , dA+ = {xedA : x x = £} , 

we would naturally be interested in the spanning clusters, which are the 
(maximal) connected clusters linking 9A__ with dA + . 

While the changes on the microscopic scale are gradual, at the perco- 
lation threshold a drastic transition is observed on the macroscopic scale, 
where the following is seen with probability which tends to 1 as a — > (for 
suitably chosen constants Const.): 
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p < p c => there are no spanning clusters; the diameters of the connected 
clusters do not exceed Const. a\ loga|. 

P > Pc => there is a unique spanning cluster, which covers the region 

"densely": its spherical voids are all smaller than Const. a\ loga|. 

At the critical point we find: 



i. the spanning probability does not vanish, at least for I < 1/2 (the 
restriction is not needed for 2D, otherwise it reflects just a limita- 
tion of the existent proof), 

ii. for each < s < 1 - 



where the bound is uniform in a (!), 

iii. there is positive probability for more than one spanning clus- 
ter in A. 

Remarks: The proofs of the above assertions are based on a number of 
different results, some of which were presented most completely for the lat- 
tice - rather than continuum, models. The behavior at p < p c follows from 
the exponential decay of the two point function in the subcritical regime 
([25,26]), the bound on the voids in the percolating phase uses the coinci- 
dence of the critical point with the limit of the slab/quadrant percolation 
thresholds ([27,28]), as explained in [12]. The statements concerning p = p c 
in arbitrary dimensions are proven in [12]. The picture described there fits 
well in the axiomatic description of the critical regime proposed in ref . [29] . 
The non-uniqueness of the spanning clusters in the two-dimensional case 
has been familiar to those following the rigorous arguments since the work 
of Russo [30] and Seymour and Welsh [30,31], but it became appreciated 
in some of the physics community only rather recently [11,13]. 

Many (though not all ) of the features of the scaling limits of critical 
percolation are expected to be shared by systems which differ on the mi- 
croscopic scale. However, beyond the general description presented above, 
the behavior at p = p c exhibits interesting dimension dependence. 



The effective conductance is an example of a quantity whose critical power law 
shows difference between some of the models listed above [32] 



P = Pc 



Prob. 




(4.2) 
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a) low dim. 



b) high dim. 



Fig. 3. Two different types of critical behavior (schematically): I) In 2D one sees 
only a finite number of clusters of the size of the volume, typically none of the clusters 
intersecting a given cube reaches far beyond it - on the corresponding scale, and in 
the scaling limit the spanning clusters are nowhere dense. II) In high dimensions (d > 
6 ?) spanning clusters proliferate (0(L d_6 )) and extend increasing distances. While 
individual spanning clusters behave as D = 4 dimensional trees, their union becomes 
dense in the scaling limit. 

5. Type I and type II critical models. The critical behavior is 
better understood in dimensions d — 2, and d > d u . c . where e? u . c . is the 
upper-critical dimension whose value for percolation is apparently d u . c , = 6 
([33-37]). It was, however, recently realized that two cases differ in ways 
which have pronounced effects on the basic structure of the continuum 
limits. 

In d = 2 dimensions, the "Russo-Seymour- Welsh theory" [30,31] (con- 
veniently summarized in [38]) can be used to deduce that: 

1. A typical configuration has only finitely many spanning clusters - 
in the sense that 

Prob Pa , a (there are in A n distinct spanning clusters) < K{n) , 



uniformly in < a < 1, with K(n) — ► (for n — > oo). More explicitly, 
K(n) was recently shown to be bounded above (C) and below (C) by 



2. Any given internal site in A is typically surrounded on an infinite 
sequence of scales (2 - ™; n = ki,k 2 , . . .) by "dual circuits" which 
separate it from all the spanning clusters, and in the limit a — > 
from all clusters whose size remains visible on the macroscopic 
scale. Of course, exceptions to the rule are found along the span- 
ning clusters, which form random "fractal" sets of lower dimension 



(5.1) 



Ae -C[C] 



[12]. 
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(< d= 2). 

In contrast, for dimensions d > 6 under an additional assumption, that 
i] = (as predicted by physical theory [34] which has been supported by 
the partial rigorous results [36,39]), we find the following behavior [12]: 

1. the number of spanning clusters grows, typically, as a~( d ~^ (as 
predicted in ref. [40]) 

2. the spanning probability tends to 1 and, furthermore, for any fixed 
open set B C A the probability that a spanning cluster intersects 
B tends to 1 

3. the diameter of the maximal cluster intersecting a given open set, 
of fixed size on the macroscopic scale ([— s, s] d ), diverges as a — > 0, 
typically being at least as large as s( d_4 )/ 2 • a~^ d ~ Q ~°^/ 2 . 

One could add that in the latter case the clusters have predominantly tree 
characteristics, and behave as D = 4 dimensional objects, as was first 
suggested in ref. [41,42]. 

The above examples, and the results presented below, motivate the 
distinction which was made in ref. [12] between the following two types of 
critical behavior. 
Type I models: The function 

„s ,. „ , ( the set [— s,s] d is connected \ 

(5.2) hmsup Prob^ a ^ ^ ^ q{ 1]d j 

.. „ , / the set [— 1, l] d is connected \ 

limsup Frob Pc>a y tQ bQundary Qf j_ 1/S; 1/g]d J 

is strictly less than one, for some s < 1. This implies lim^oo h(s) = 0, 
which means that there is no percolation in the scaling limit. 
Type II models: 

, . / the set [— s, s] d is connected \ 

(5.3) Frob Pc . a y tQ boundary of yd J — » 1, 

for any < s < 1. 

Remarks: 1) Presumably all other behavior is ruled out for the models 
considered here, but that was not proven. 

2) The function h(s) is obviously submultiplicative, i.e., satisfies h(si-s 2 ) < 
h(si)h(s 2 ) for < si,s 2 < 1- Standard arguments imply the existence of 
the limit, and its positivity for Type I systems: 

(5.4) lim^M = A>0. 

s\0 logs 

In other words, in Type I critical models h(s) = s x+ °^ for s — > 0, with 
some A > 0. In terms of the standard, though not yet fully proven picture 
[4], the exponent A is related to the Hausdorff dimension of the Incipient 



~h(s) 
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Spanning Cluster (or of the ISC), d w , as: 

(5.5) A = d — d w . 

In this article we shall not discuss the scaling limit of models of of Type 
II (i.e., the case d > d u . c .). Not that this would be uninteresting: Hara and 
Slade proposed (as mentioned in [43]) that the limit for individual clusters 
is related with the Integrated SuperBrownian Excursion process (ISE) of 
Aldous [44]. Furthermore, looking at the ensemble of all the macroscopic 
scale cluster we find that in a sense which still has to be made explicit (the 
one presently in mind is a very weak one) one may anticipate two surprising 
features: percolation at the critical point, and infinity of infinite clusters 
[12]. 

Remark: The above observation may alarm those familiar with lattice 
percolation models, but the apparent conflict with the general uniqueness 
Theorems [45,46] is not a real contradiction; the general result of Burton- 
Keane [46] requires discreteness on some short-scale. 

Our discussion of the scaling limit continues now in the generality of 
of Type I models. This covers the case of 2D, and presumably applies also 
to dimensions d = 3, 4, 5 - though there are no rigorous results to support 
such claim. We of course limit now our attention to the critical regime. 

6. Formulation of the scaling limit; type I models. 

6.1. The question. The quantitative description of the continuum 
limit is expressed through a number of functions. Two examples are de- 
scribed in Figure 4. The pictures drawn there refer to events defined on 
the macroscopic scale, with the probabilities considered in the limit a — > 0. 

The question we shall address next is what mathematical object would 
capture, in a natural way, those geometric features of the Incipient Span- 
ning Clusters which are visible in the scaling limit. This can be rephrased 
as asking what stochastic-geometric object embedded in R d can be associ- 
ated with the functions referred to above, and others of this kind, so that 
they can be naturally viewed as the connectivity probabilities (and not just 
limits of . . . ). 

In the absence of a direct insight a canonical approach could be to to 
define the object by the list of its quantifiers, and introduce on the space 
of those some minimal a - algebra which would allow to bring in proba- 
bilistic notions. However, as the example of the Wiener process (Brownian 
motion) shows, it may be worthwhile to learn the regularity properties of 
the geometric object under consideration. 

6.2. The first attempt: ISC as a subset of A. 

At first sight, it seems natural to regard the ISC as closed subsets of 
the region A C R d . However, this formulation will not serve our purpose. 
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G (x , x , x ) = Prob 

D, e 1 2 ' 3 J 




Fig. 4. Two functions associated with the continuum limit (a = (lattice spacing) — ► 0) 
of Type I critical models: 1) F is the limit of the probabilities that boundary segments 
are connected as indicated, 2) G is the limit of probabilities that the neighborhoods of 
the points {x\ , X2, xz} are connected. 

Before dismissing this attempt, let us recall that what makes the col- 
lection of closed subsets of A C R d into a particularly convenient space for 
the description of random fractals ([47]) is the following classical result. 

Theorem 1. For any compact metric space A, the space 



is compact in the corresponding Hausdorff metric. 

The Hausdorff metric h(A, B) is defined so that h(A, B) < e if and 
only ifdist(x, B) < e for any x G A, and dist(y, A) < e for any y E B. 

Alas, the features which are of interest to us, such as the existence and 
location of connecting paths, are not continuous in the Hausdorff metric. 
Moreover, the configurations of critical percolation models are typically 
among the points of discontinuity. The reason is that in typical configu- 
rations there are "choke points" where a small-scale change, possibly of a 
single bond, or cell, would drastically alter the available connecting paths 
(as indicated in Figure 3a, and more explicitly in Figure Id of ref. [12]). 
Such a change shifts the point in 7i(A) only a distance of the order O(a), 
which is not detectable in the scaling limit. However, the effect on the 
available connecting paths is clearly visible on the large scale. 

6.3. Holder continuity of the connecting paths. Since the scal- 
ing limit of the set does not capture the information on the realized connec- 
tions, it is natural to include that information explicitly in the description 



(6.1) 



H(A) = {A c A : A is closed } 
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of limit. Some of it is expressed through the collection of the realized self- 
avoiding paths, each given by a continuous function g : [0, 1] — > A. In the 
terminology of [4], we are including both the backbone (BB) paths and the 
paths connecting the dangling ends to the backbone. 

Potential obstacles in describing the configuration through the realized 
paths are: 

1. The possibility that as the short scale is refined the connecting 
paths could become increasingly irregular. It is not a-priori obvious 
that in the limit a — > the connections can still be expressed 
through continuous functions. (One could worry here about the 
need to consider more general continua [connected closed sets]. 
Their collection is somewhat unwieldy, e.g., some continua do not 
support the image of any continuous non-constant function.) 

2. It is not initially clear whether the information provided by the set 
of the connected paths suffices for questions concerning higher or- 
der connections. If not, then one might need to list also connected 
line graphs of higher complexity. 

The first concern is completely answered by the following result [48,49] 
(see Note Added in Proof, next page). 

Theorem 2. For any critical Type I percolation model all the realized 
connected (self- avoiding) paths in a compact region A C R n can be simul- 
taneously parameterized by uniformly continuous functions, g{t) < t < 1, 
satisfying the Holder continuity condition: 

\g(ti) -g(t 2 )\ 

(6.2) ~ / < k(u) , for all < t u t 2 < 1 , 

with some fixed < A < 1/d and a configuration dependent continuity 
modulus for which 

Prob Pc , a {k(uj) >t) < g(t) 

(6.3) -> , for t -> 00 

uniformly in a. 

In other words: in Type I models, in the critical regime one seldom 
finds a connected path in A which cannot be "traced" in a unit of time by 
means of a "fairly regular" function. The continuity condition we use is 
consistent even with a "fractal" landscape, and consequently the regularity 
does not deteriorate as a — ► 0. (That is not true for p > p c .) 

The self-avoidance condition is applied only on the microscopic scale, 
and should be interpreted in the sense which is natural for the model, e.g., 
for bond percolation the paths should not repeat any bond, and for the 
random Voronoi-tessellation the paths should not re-enter any cell. We 
note that the paths need not appear self-avoiding when viewed on the 
macroscopic scale (in the limit a — > 0). 
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Concerning the second of the above reservations, the question has a 
simple answer in d = 2 dimensions, though the situation in higher d is still 
not as clearly resolved. The basic issue is: how to determine if a pair of 
paths which on the scale of the continuum seem to intersect are actually 
connected on the microscopic scale. As we discussed, there are situations 
in which two connected paths come within distance a without touching 
and without there being another path linking the two. Conveniently, at 
least in d = 2 dimensions such close encounters 1st kind do not occur at 
non-terminal points, in the sense which is stated precisely in ref. [12], and 
proven in ref. [48]. 

6.4. Our choice: ISC as a collection of realized Holder 
continuous paths. 

To formulate the limiting representation of the Backbone and the In- 
cipient Spanning Clusters, we find it convenient to first represent the perco- 
lation configuration by the random collection of all the realized (connected) 
paths in A which are regular and self-avoiding in the sense explained above. 
That random object, which we call the percolation Web (W, or W(co)), is 
of the form 

(6.4) W cC 1/d (lO,l],R d ) , 

where C 1/d ([0, 1], R d ) is the space of Holder-continuous functions: 

C 1/d ([0,l], J R d ) = (. 9 GC([0,l], J R d ) : sup ^ ~ f f < oo) 

I 0<ti,t 2 <l 1*1 - *2| ' J 

(6-5) 

We shall also denote by W\(u)) the restriction of W{w) to functions with 
range in A C R d (i.e. to C([0, 1], A)). 

In d = 2 dimensions, the range of values of W{u) is constrained by 
the following consistency conditions, of which the first three are obvious, 
but the fourth one reflects a non-trivial fact (there [typically] are no close 
encounters of the 1st kind, in the limit a — > 0). 

Percolation web consistency conditions 

CI (Closure) W\(u) is closed as a subset of C([0, 1], A). 

C2 (Reparametrization-invariance) For each realized path / £ Wa(u>) 

any path of the form f(t) = /(r(i)), with r(-) continuous and of 

bounded derivative [or just a Lipschitz function], is also realized. 
C3 (Splicing stability) If two paths of W\(lu) intersect at non-terminal 

points, then the paths obtained by different "splicings" of the four 

resulting segments are also realized. 

Note Added in Proof: An equivalent but possibly more appealing formulation 
is to describe ISC through the collection of the realized paths of finite tortuosity. This 
approach, along with some basic results concerning random curves with bounded tortu- 
osity, is being developed in a joint work with Almut Burchard [49]. 
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We denote by f2 A the collection of subsets of C([0, 1], A) satisfying the 
consistency conditions CI - C3. 

(Trying not to be too formalistic here, let us just note that f2 = fl R d 
is a complete separable metric space. The continuum limit of critical per- 
colation models in a macroscopic region A will be described by probability 
measures on the natural a algebra on 0(A).) 

The Backbone can now be described by the collections of paths in the 
Web which traverse A, 



and the collection of the Incipient Spanning Clusters is described by a 
collection of pairs of paths, of the form: 



where (/, g) £ S\(uj) is taken to imply that there is an actual contact 
between the paths (at the microscopic level, which is otherwise no longer 
visible) . 

7. The web — an existence result. From the perspective of the 
continuum limit, the microscopic model is a construction scaffold. When it 
is removed, a more remarkable structure is exposed (as in Emily Dickinson's 
metaphor) . 

Once we have a direct way to formulate the continuum theory, it is 
mathematically natural to restart the discussion and pose the question in 
the standard existence and uniqueness terms. It should be appreciated that 
a successful formulation of the uniqueness question will shed light on the 
phenomenon of universality in critical behavior. 

Following is an existence result [48] . 

Theorem 3. For each dimension in which the critical behavior is Type 
I, there is a one-parameter family of probability measures (fi) on fi which 
have the following properties. 

1. (Independence) For disjoint closed regions, A C B C R d , Wa(w) 
and Wb(^) are independent [as random variables]. 

2. (Euclidean invariancc) The probability measure is invariant under 
translations and rotations. 

3. (Regularity) The spanning probabilities of compact rectangular re- 
gions are neither nor 1 : 



(6.6) 



BA(w) = {/eW A (w):ffi(0) = - 



(6.7) S A (u) = Uf,g) 



f 6 Wa(u), jG B(u>), with 
1(0) = g(t) for some and < t < 1 



(7.1) 




there is in [—s,s] d a 
left right spanning cluster 



> 
< 1 



A convenient parametrization within the family of measures is the 
crossing probability R\ . 
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The measures are constructed as continuum limits of sequences of mod- 
els with suitably adjusted percolation densities. If the standard picture is 
correct, the density needs be adjusted with the lattice spacing as: 

(7.2) p = p c + (R 1 - 1/2) Const. a 1/v . 

In order to end with a rotation invariant measure, we start from either the 
droplet percolation model, or the Voronoi-tessellation percolation model. 

Some open problems: 

Convergence — A characteristic shortcoming of the available methods is 
the lack of proof of convergence of the scaling limit. Our con- 
struction relies on compactness arguments, which guarantee con- 
vergence along subsequences. Proof of convergence will be an 
outstanding technical contribution to the subject. Alternatively 
stated, this is a question of uniqueness of the scaling limit. 
Uniqueness — A broader formulation of the uniqueness question is: 

do the three conditions seen in the existence result, Theorem 3, 
limit the range to only the one-parameter family of measures? 
If not, are there additional assumptions which would? 
Also: is the full rotation invariance of such measures implied by 
just the rotation invariance of the cube-crossing probabilities? 
Since the measures in question include all the continuum limits of critical 
percolation models, positive answer would cast in a clear mathematical 
form some of the expected universality of critical behavior - in a sense 
which was clearly articulated only relatively recently, in Langlands et.al. 
[5] and related works [15, 6-9]. 

Other mathematical challenges arc to affirm (or test) the Rcnormal- 
ization Group picture, which suggests some exact properties for the con- 
structed measures, and establish the conjectured conformal invariance of 
the critical measure (a special member of the one-parameter family). 

8. Relation with the renormalization group. The continuum ob- 
ject constructed in Theorem 3 bears an interesting relation to the RG 
picture. While still no sensible formalism has been found for an exact 
representation of the renormalization group transformation as a map, the 
one-parameter family of measures presented in Theorem 3 may be viewed 
as corresponding exactly to what would be the unstable (i.e., expanding) 
fiber extending from the critical point in such a space. Along this fiber, 
the RG maps coincide with dilatations. 

It is interesting to plot the joint values of R(s) and R(l) (say, for s = 2). 
The standard RG picture leads us to expect there the S*-shape of a function 
which as a map of the unit interval into itself has one unstable fixed point 
and two stable fixed points at and 1 (as in Figure 5). The slope of the 
function at the unstable fixed point should be exactly s 1 ^ . (For 2D the 
predicted value is v = 4/3, den Nijs [50]). 
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Fig. 5. The conjectured relation between two crossing probabilities in the one- parameter 
family of the constructed measures (schematically). The map this function generates on 
[0, 1] corresponds to the action of the Renormalization Group along the expanding fiber. 
Assuming the standard RG picture, the slope of the function at its middle [unstable] 
fixed point is exactly 2 1 /", v being the correlation-length critical exponent. 

Even in the absence of an exact setup, the renormalization group pic- 
ture has provided very effective approximate tools [51]. One may view 
Figure 5 as presenting a limit of the "cell to cell" renormalization group 
map, which recently became appreciated as one of the more effective ap- 
proximate RG methods [6,52,53]. 

9. Fractal structure. In Type I models, the number of clusters con- 
necting the boundary of a cube A c R d 1 centered at 0, with the contracted 
cube sA, < s < 1, is finite; in the sense that the probability for observing k 
such clusters obeys bounds which are uniform in the short-distance/lattice- 
spacing a (as a — > 0), and decay to zero for k — ► oo. Correspondingly, the 
scaling limit exhibits (a.s.) only a finite number of such clusters (once one 
knows how to count them), and altogether only countably many macro- 
scopic size clusters in A (the infinity is caused by the union over all scales) . 
Furthermore, the number of "left <-> right" spanning clusters remains finite 
in the continuum limit. 

Details of the expected fractal structure were discussed in the review 
article of Stanley [4] in the context of lattice models. New considerations 
are added when one looks at the scaling limit. To present some basic 
results, let us start with- 

Definition: 1. For a given configuration of the percolation Web, W\(u>), 
we say that two sites, x,y € A are connected if W\(u>) includes a path 
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visiting each of them. 

2. The connected cluster of a site x G A is the union of the sites y G A 
connected to it. We denote it C(x) = C u {x). 

3. C(dA) is the collection of sites connected (as in 1.) to the boundary dA. 
4- The ramification number R u (x) is the maximal value of k for which 
there are k paths in W(u>) starting at x and otherwise non-intersecting. 

There are some surprises: first is the lack of transitivity of the relation 
"x is connected to y" . We view this not as a shortcoming of the terminology, 
but rather as an expression of an interesting phenomenon, related to the 
existence of tenuous connections (when a pivotal bond is reversed, one 
obtains a configuration in which two realized paths meet without being 
connected on the microscopic level, an example is indicated in Figure 3a). 

The second "surprise" is the first statement in the following list of 
properties of the scaling limit, ref. [48]. 

Theorem 4. For Type I critical models, in typical realizations of the 
percolation Web: 

1. The connected clusters of (Lebesgue-)ahoaost all sites x G A contain 
no other site, i.e. R{x) = or 

(9.1) C{x) = {x} . 

2. The collection of sites violating eq. (9.1), which includes the random 
set C u {d\), is of Hausdorff dimension < (d — A) (as defined by eq. (5.4). 
Furthermore, 

3. The above set is of finite ramification; there is a non-random value 
k < oo such that R hJ {x) < k for all (not just a.e.) x G A. 

For 2D, we guess that the maximal ramification number is about k = 5, 
though that is still not fully resolved. 

10. Conjectured conformal invariance. The particular measure 
which corresponds to the fixed-point value of Ri is expected to be fully 
dilatation invariant. It is conjectured that it is also strictly covariant under 
conformal maps. 

The conformal invariance is expected in any dimension, but the con- 
formal group is particularly large in 2D and it is there that the related 
considerations have been shown to have powerful consequences, and have 
led to explicit predictions concerning critical behavior in a rich collection 
of models [54,55]. For the rest of this section we restrict the attention to 
two dimensional systems. 

The 2D scaling limit was studied numerically by Langlands et.al. 
(LPPS) [5]. In addition to testing the universality of the spanning prob- 
abilities (in a sense which broke new grounds [6,7,8,9]), LPPS asked how 
should the aspect ratio of a parallelogram be adjusted with the angle (see 
Figure 6), if one wants to keep the spanning probability constant. This au- 
thor's suggestion that the criterion should be conformal equivalence fitted 
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Fig. 6. The setup for the question of Langlands et al. [5] and Cardy's surmise. 

well with the numerical data. A much more complete answer was proposed 
by J. Cardy [15], who produced a differential equation for the upper-half- 
plane version of the problem (to which the original one is reduced via the 
corresponding Riemann map) . Cardy's equation for the crossing probabil- 
ity has a unique solution with the natural boundary conditions, and the 
resulting function of the angle a and the length L (see Figure 6) was found 
to be in perfect agreement with the numerical results of LPPS. 

Cardy's equation drew on a field-theoretic perspective, and on an ex- 
trapolation based on relations with some other models of Statistical Me- 
chanics. Other insights followed, including Pinson's proposal [56] for the 
exact values of the probabilities of different windings for twisted boundary 
conditions. The reader is referred to [16] for an account of these develop- 
ments. 

The bottom-up derivation of the conformal invariance is still an open 
challenge. There are some partial results, such as the following one (whose 
proof is not that difficult) which is conditioned on a strong assumption 
about the scaling limit. 

The result stated below follows [48], but we note that a related state- 
ment with a somewhat different formulation was presented in [23] ([the 
two were arrived at independently]). To formulate the proven assertion, 
let us first present a statement which upon some consideration appears 
believable, though its proof has still eluded us. 

Conjecture: The Voronoi-tessellation percolation models with position- 
dependent density profiles of the form p(x) = t ■ g(x), with g(-) continuous 
and non-vanishing in A C R 2 , the limit t — > oo exist (for the ISC and for 
the Web processes), and is independent of the function g(-). 

Theorem 5. In d = 2 dimensions, under the above conjecture and 
assuming also the critical behavior is Type I, the scaling limit of the Voronoi 
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percolation model is conformally invariant, in the sense that for any map 
T : A — > R 2 which is invertible and conformal on A, the image of the ISC 
[Web] process in A coincides with the ISC [Web] process in TA. 

A particular expression of the conformal invariance is that the function 
defined in Figure 4 satisfies 

(10.1) F D ({ 7l ,..., 7fe }) = F TD ({ 7l ,..., 7fe }) , 

where {71, . . . ,-fk} refers to a collection of boundary segments which are 
to be connected by a common cluster. (One may of course generate a large 
number of other, similar, functions.) 

One may employ also the Riemann map which takes the interior of D 
conformally onto the upper-half-plane; Figure 6 indicates two such events 
for which conformal invariance implies equality of probability. In case the 
boundary of D consists of four segments, conformal invariance implies that 
the probability depends only on one cross-ratio, i.e., is given by a function 
of the form 

(10.2) f(xi, ...,Xi) = 4>{ u(xi, ...,x 4 )), 
with 

(xi - x 2 )(x 3 - Xi) 



(10.3) u(xi, . . . ,Xi) 



{xi - x 3 )(x 2 - x A ) 



11. Cardy's equation. The equation which Cardy proposed for the 
above quantity can be transcribed as [15,16]: 



with the boundary values </>(0) = 0, <j)(l) = 1. The solution can be presented 
in the integral form: 

(H-2) 0(u) = f ' — I Norm. 



[x(l-x)} 2 / 3 

Some algebraic aspects of Cardy's equation are more visible when it is 
expressed in terms of the differential operators: 



+ : - 

Ox, 



(11.3) £ n = -J2( X j~ X ^ 

In terms of these, the equation (as originally presented) is 

(11.4) {C. 2 -lC 2 _ l )f{x u ...,x i ) = . 

Cardy's original argument [15] invoked an analogy with equations de- 
scribing the effects on the free energy of certain Potts models of changes 
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in the boundary conditions. Various step in this approach are still beyond 
the reach of rigorous methods. Let us, however, point out that the very 
intuitive hypothesis of invariance under conformal maps which preserve the 
relevant domain (plus the necessary differentiability) implies: 

(11.5) A (£-2 - f(x u ■ ■ ■ , a*) = . 

The derivation of this equation is an elementary and amusing exercise in 
the Virasoro algebra, of the commutation relations: 

(11-6) [C n ,Cm] = {n-m)£ n+m , 

starting from the observation that the conformal invariance assumption 
implies 

(11.7) £>/(•■■) = ,and - . 

Equation (11.7) holds since the two operators generate flows preserving 
/(...), being associated with Mobius transformations which preserve the 
upper-half-plane and the point X\ (£-i, however, corresponds to a shift 
of x 2 , ■ ■ ■ which leaves Xi behind). Thus, when the two factors in eq. (11.5) 
arc transposed /(. . .) is annihilated by C\. The commutator of the two 
terms is a combination of Co, C\ and a term proportional to C-i which is 
eliminated through the judicious choice of the coefficient 3/2 in eq. (11.5). 

The transition from eq. (11.5) to eq. (11.4) amounts to the removal of 
the former's leftmost factor. That is however a big step, since it transforms 
an equation with limited content into one which completely determines the 
solution. (Although one should not be too dismissive of eq. (11.5): in 
terms of the conformal invariance structures, it conveys the fact that for 
the percolation problem the conformal charge is zero, which is a speculation 
born out by the Russo-Seymour-Welsh theory [30,31].) 

Though we are still short of a proof, it turns out that one may explain 
from the Continuum Percolation Web's perspective a mechanism for the 
extraction of Cardy's equation, eq. (11.4), from eq. (11.5). The argument 
employs a plausible description of an effect caused by the separation of 
scales ([48]). 

12. Selected characteristic exponents. An interesting aspect of 
the explicit solution, cq. (11.2), is that it yields the power e 1 /' 3 for the 
probability that a short interval, of size e, is connected in the half plane a 
distance of order 1 away. This exponent does not appear to be obvious to 
the naked eye. In Figure 7 we present it next to two other characteristic 
exponents which have a more elementary derivation (discussed along with 
other examples and applications in ref. [21,10,48]). 

Less one would be lulled by the simplicity of the exponents seen in 
Figure 7, let us mention that the probability of a full annulus to be spanned 
is predicted to behave as (r/i?) 5 / 48 , den Nijs [50]. 
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Fig. 7. Some characteristic exponents associated with spanning probabilities in 2D. 
Solid lines indicate spanning clusters, and dotted lines indicate spanning dual, separat- 
ing, clusters. Example a) is implied by Cardy's equation, b) and c) are based on more 
elementary observations, [21,48]. 

Question: What are the values of the characteristic exponents for k dis- 
joint spanning clusters, in the full or cut annulus, for all the other values 
of fc? 

The exact answer may be within the reach of the methods of Conformal 
Field Theory [55,15], or the Coulomb-gas representation, [57], which has 
been recently discussed by T. Spencer[58]. It was recently proven by other 
methods that for large k the exponent is of the order of Const, k 2 , [12]. 

13. Relation with field theory. The continuum object described 
here in earlier sections is related to a number of "field theories" . Cardy's 
equation for the function F drew on its analogy with the vacuum expecta- 
tion value of a product of operators switching the boundary conditions (in 
Potts models). Another field theoretical object is related with the function 
G, also defined in Figure 4. Its simplest manifestation is in the one-point 
function defined next. 

In the scaling limit of a Type I model, the probability that a given 
site x <E D is connected to the boundary of the domain D is zero, but 
the probability that a small ball around it, B(x,e) — {y : \y — x\ < e}, 
is connected to 3D is positive. Based on the ideas mentioned in Section 5, 
one expects that for a suitable d w the following limit exists 

(13.1) h D (x) = lim e - fe ( 2 - d ») Prob(B(x,e) is connected to 3D) . 

If so, and if the probability measure maps covariantly under conformal 
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maps, then, for any conformal transformation T : D — ► R d 

(13.2) h D {x) = {\&et{dT/dx)\ 1 l d f- d ™ h TD {Tx) . 

The function h d {x) can be viewed as the expectation value of an entity 
defined by the limit 

(13.3) ip D (z) = w- lim e -( 2 - d ») l[B(x,e) is connected to 3D] . 

The limit is initially interpreted in the weak sense (i.e., expectation values 
of . . . ). If we try to give ^d(z) a stronger meaning, in the context of the 
continuum process discussed in earlier sections, we find that for a given 
realization of the continuum percolation Web, iPd( z ) is zero at almost all 
points, {z G D}, but it diverges along the fractal set connected to the 
boundary. It would be natural to think of it as a distribution-valued ran- 
dom field, but clearly some more thought should be given for a complete 
development of this interpretation. Related concepts on the horizon arc: 
the stress-energy tensor for percolation, and operator product expansions 
for entities like *Pd{z) (discussed in other contexts in [59,54,60,55]). 

Let us conclude by noting that similar considerations apply to the n- 
point function associated with the function Go,e defined in Figure 4. 

If, as in eq. (13.1), the following limit exists, 

(13.4) g D (zi,...,z n ) = lime- k ( 2 - d ^G D , e (z u ...,z n ) , 

then under the "plain" conformal invariance hypothesis the resulting func- 
tion should satisfy 

n 

(13.5) g D (zi,...,z n ) = H\T'( Zj )\ 2 - d - g TD (T Zl ,...,Tz n ) . 

We switched here to notation appropriate for 2D, where a rich class of 
conformal transformations is provided by the maps T : D — > C which are 
analytic and invertible on D. Equation (13.5) is reminiscent of the transfor- 
mation law of the vacuum expectation values of products of field operators 
in conformal field theory [54] , although the higher order connectivity func- 
tion were not yet transcribed into such expectations of products of local 
field operators, and percolation seems to be missing the reflection positiv- 
ity (corresponding to the positivity of the inner product) which plays an 
important role in field theories. 

It would be interesting to see further development of the field theoretic 
content of the continuum theory described here. This should be done both 
within the perspective native to percolation models, and in the direction 
of links with other field theoretic structures. Existence of such links is 
suggested by the Fortuin-Kasteleyn relations of percolation with Ising and 
Potts models [61], and the relation of the latter with (j) 4 and other field 
theories. 
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